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ABSTRACT 

It is generally believed that the radiation of relativistic particles in a curved magnetic field proceeds 
in either the synchrotron or the curvature radiation modes. In this paper we show that in strong 
curved magnetic fields a significant fraction of the energy of relativistic electrons can be radiated 
away in the intermediate, the so-called synchro-curvature regime. Because of the persistent change 
of the trajectory curvature, the radiation varies with the frequency of particle gyration. While this 
effect can be ignored in the synchrotron and curvature regimes, the variability plays a key role in 
the formation of the synchro-curvature radiation. Using the Hamiltonian formalism, we find that 
the particle trajectory has the form of a helix wound around the drift trajectory. This allows us to 
calculate analytically the intensity and energy distribution of prompt radiation in the general case 
of magnetic bremsstrahlung in the curved magnetic field. We show that the transition to the limit 
of the synchrotron and curvature radiation regimes is determined by the relation between the drift 
velocity and the component of the particle velocity perpendicular to the drift trajectory. The detailed 
numerical calculations, which take into account the energy losses of particles, confirm the principal 
conclusions based on the simplified analytical treatment of the problem, and allow us to analyze 
quantitatively the transition between different radiation regimes for a broad range of initial pitch 
angles. These calculations demonstrate that even very small pitch angles may lead to significant 
deviations from the spectrum of the standard curvature radiation when it is formally assumed that 
a charged particle moves strictly along the magnetic line. We argue that in the case of realization 
of specific configurations of the electric and magnetic fields, the gamma-ray emission of the pulsar 
magnetospheres can be dominated by the component radiated in the synchro-curvature regime. 
Subject headings: gamma rays: general - magnetic fields - radiation mechanisms: non-thermal 


1. INTRODUCTION 

The energy distribution and intensity of radiation of 
a charged particle moving in the magnetic field with a 
relativistic speed is determined by the curvature of its 
trajectory. The trajectory in the curved magnetic field 
has a structure of a helix wound along the so-called drift 
trajectory , which describes the particle motion averaged 
over the fast oscillations. Thus, the curvature of the tra¬ 
jectory is formed by the curvature of the helix itself and 
the curvature of the drift trajectory; they are determined 
by the strength and the curvature of the magnetic field, 
respectively. If the main contributor is the curvature of 
the helix (i.e., the bending of the drift trajectory can be 
neglected), the radiation is defined by the strength of the 
magnetic field. This is the case of the synchrotron radi¬ 
ation. On the other hand, if the helix is stretched (i.e., 
the particle moves along the drift trajectory with a little 
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wobbling around it), the main contribution to the radia- 
tion is due to the curvature of the drift trajectory. This is 
the case of the curvature radiation, which is determined 
by curvature of the drift trajectoryQ. Obviously, the syn¬ 
chrotron and curvature radiation regimes are approxima¬ 
tions that work with a high precision as long as the cur¬ 
vatures of the helix and the drift trajectory differ signif¬ 
icantly. However, when they are comparable their joint 
action may result in features that are significantly dif¬ 
ferent from those described by the standard synchrotron 
and curvature radiation formalism. These two contribu¬ 
tions to the curvature can reinforce or extinguish each 
other, and therefore cause a significant variability of the 
curvature of the particle trajectory. It results, in turn, 

1 Note that usually in the literature the curvature radiation is 
considered as radiation when particles move along the magnetic 
line. Strictly speaking, this is a formal but, in fact, not a correct 
assumption because a particle can move only along a drift trajec¬ 
tory. However, because the drift trajectory is close to the magnetic 
line in the strong magnetic field, this assumption does not lead to 
a significant inaccuracy. 
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in a radiation that varies with the frequency of the par¬ 
ticle gyration. The observer detects the time-averaged 
energy distribution of this promptly variable radiation, 
which appears rather different from the spectra of emis¬ 
sion produced in the synchrotron or curvature radiation 
regimes. 

A few attempts have been made in the past to study 
thi s intermediate re g ime o f radiation, which we, follow¬ 
ing [Cheng^2Ein3 (Il996h . call synchro-curvature radi¬ 
ation. However the approaches and assump tions used in 
the p r evious treatments of the problem (|Cheng fc Zhand 
fl99?l IHarko fc Chmel l200l are generally not se lf- 
consis tent. For example, in the paper of lCheng fe Zhanal 
(jl996t ). where the particle motion is introduced a priori 
(by hand), the square of t he velocity is not conse rved. 
The more recent results of IHarko fc Cheng! (12 002 1 also 
seem not correct. In particular, the acceleration in this 
work is not variable as it follows from the solution of 
equations that describe the particle trajectories (com¬ 
pare the formula given by Equation (5) in their paper 
and the expression given by Equation (1401) in the present 
work). On the other hand, we cannot fail to mention the 
apparently underrated work of lSobolev ( 2006 1, where the 
author has thoroughly treated the problem of radiation 
of particles in the magnetic field. Some principal results 
and conclusions of this paper are in a good agreement 
with the results of Section [2] of the present work, al¬ 
though they have been obtained using a rather different 
computational approach to the problem. However, these 
results are of rather theoretical interest because the en¬ 
ergy losses of particles are ignored. 

In this work, using the Hamiltonian formalism, we 
solve the equations of the motion of a charged particle in 
the azimuthally symmetric magnetic field. This allows 
us to self-consistently describe the local properties of the 
particle trajectory with a time-dependent curvature, and 
calculate the time-averaged spectrum of radiation in an 
arbitrary magnetic field. The analytical solutions de¬ 
rived in an approximation that ignores energy losses of 
electrons provide a useful tool for understanding of dis¬ 
tinct feature of the synchro-curvature radiation, and an¬ 
alyzing the conditions of its transition to the limits of 
synchrotron and curvature regimes. 

On the other hand, for the quantitative description of 
propagation and radiation of particles, one has to treat 
the energy losses of electrons properly. For this reason, 
we have performed numerical integration of the equa¬ 
tions of motion in the dipole magnetic field, taking into 
account the radiation reaction force. The dipole mag¬ 
netic field can be used as a good approximation for the 
strong magnetic held in the vicinity of compact objects 
like neutron stars. The results of numerical calculations 
allow us to study the evolution of radiation for differ¬ 
ent initial pitch-angles. These results demonstrate that 
even a small initial deflection of the motion of electrons 
from the drift trajectory could lead to dramatic devia¬ 
tions from the radiation features expected in the con¬ 
ventional curvature radiation regime. The question of 
characteristic values of pitch-angles in the pulsar mag¬ 
netospheres remains an open issue. The answer can be 
obtained if one brings the electron acceleration into con¬ 
sideration. However this requires a knowledge of the rel¬ 
ative arrangement of the electric and magnetic fields. We 
will discuss this issue in the context of the acceleration 


of electrons in the simple model of the electric field. 

This paper has the following structure. In Section [2] 
we analytically solve the equations of the motion of a 
particle in the magnetic field of a constant curvature. 
Using the curvature of the trajectory, in Section [3] we 
derive analytical expressions for the radiation spectrum 
of synchro-curvature radiation. Subsequent sections deal 
with a numerical treatment of the particle radiation in 
some specific configurations that are of astrophysical in¬ 
terest. Sectionjll as well as AppendiceslAlandlBl describe 
relevant formalism for the particle motion and its radi¬ 
ation. In Section [5] we discuss the results of numerical 
calculations, namely the radiation spectra and regimes 
of radiation for different initial conditions. In Section [ 6 ] 
and Appendix [C] we investigate the impact of the com¬ 
bined action of electric and magnetic fields on the par¬ 
ticle acceleration and radiation. Finally, in Section 0 we 
summarize the main results and conclusions. 


2. PARTICLE TRAJECTORIES AT SMALL PITCH ANGLES 

The motion of a relativistic charged particle in a strong 
magnetic field is accompanied by a radiation that in gen¬ 
eral terms can be called magnetic bremsstrahlung. While 
moving in a homogeneous field, a charged relativistic par¬ 
ticle gradually radiates away its energy. However, the 
component of its velocity, which is parallel to the mag¬ 
netic field, Dm, remains unchanged. This statement is 
obvious for d,, = 0 , therefore in an arbitrary coordinate 
system, which moves along (or opposite) the magnetic 
held, we have Dy = const. Thus, 



pi 


P II 


= const, 


(1) 


where p,, and p x are the parallel and perpendicular com¬ 
ponents of the momentum, respectively. This equation 
allows us to find the energy of the particle after the radia¬ 
tive damping of the perpendicular component of motion. 
We denote the initial and final values of the energy, and 
the parallel and perpendicular components of the mo¬ 
mentum as E, P||, p_|_, and E ', Py, p' ± = 0, respectively. 
From Equation © it follows that 


P II =P|| 


\/p\ + m 2 c 2 


E' = E ■ 


y/pi + ro 2 c 2 


( 2 ) 


If the initial perpendicular momentum is nonrelativis- 
tic, p ± < me, then E' ~ E (i.e., the damping of the per¬ 
pendicular motion does not have an effect on the particle 
energy). However, for me <C Pj_ <C Py, we have E' <C E 
(i.e., the particle loses almost all its energy during the 
damping of the perpendicular motion). For instance, for 
a particle with the initial Lorentz factor 7 and the ini¬ 
tial pitch angle a = IOO /7 <C 1, the final Lorentz factor 
is reduced by two orders of magnitude, 7 ' = 7 /100. At 
p ± toc, the momenta py and p^ vary similarly with 

time because of the constant ratio Pj_/py = const (see 
Equation Jl])). The constancy of the ratio p x /py has a 
simple explanation. The ultrarelativistic particle emits 
photons in the direction of its momentum p within the 
angle < 1 / 7 , therefore in the case of I/7 <C P x /Py, the 
recoil momentum is directed opposite to the vector p. 
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No matter how small the pitch angle is, as long as it is 
much larger than 1 /y, the particle will lose a considerable 
amount of its energy. 

In the homogeneous field, the particle, after losing the 
perpendicular component of its momentum, continues its 
motion along the magnetic field line. The situation is 
quite different when the field is curved. If the velocity 
vector is parallel to the field line initially, the accelera¬ 
tion a| ( _ Q = 0. Correspondingly the curvature radiuifl 
R c L „ = c», thus no radiation is emitted at t = 0. In the 
curved magnetic field, a rectilinear uniform motion is not 
possible, so the acceleration and radiation of the parti¬ 
cle are unavoidable. For all trajectories located close to 
the drift trajectory, the acceleration and curvature ra¬ 
dius appear to be time-dependent. This implies that the 
radiation also should vary with time, thus for calculation 
of the measurable characteristics of radiation one should 
average the relevant distributions over time. 

We start with a simple model of particle motion that 
allows exact solutions. Let us consider the case when 
the magnetic field has the same symmetry as the field of 
an infinitely long straight wire. Then, in the cylindrical 
system of coordinates, (r, </>, z), the field only has an 
azimuthal component; it only depends on the distance to 
the axis of symmetry, r, but not on z and (f >. We denote 
by (e r , e^, e z ) a moving system of unit vectors linked to 
the point (r, </>, z). The magnetic field can be expressed 
as B = B(r) e^. The field lines are circles of a radius r, 
which means that at a particular distance from the axis 
of symmetry the magnetic field has a constant curvature. 
Particles in such a field can move along trajectories close 
to the circular field lines. 

The equation that describes the particle motion is 


m r yv = - (v x B). 
c 


(3) 


Let us find a solution for which r = tq = const, and the 
velocity 

v = + v z e z . (4) 


Note that the components v$ and v z are constants. The 
vector is uniformly rotating with an angular frequency 
ui = v^/tq , therefore the derivative de^/dt = —V 0 e r /r 
and the vector e z = const. Substituting Equation 
into m, we find 


v 4> e 

m 7 — e r = - v,B e r . 

r 0 c 

From here it follows that 

v z _ mcv 07 
V 0 eB ro 


(5) 

( 6 ) 


Writing the velocity components in the form v z = v sin a, 
V 0 = vcosa (the angle is counted from the (x,y) plane, 
i.e., —7 t/2 < a < 7 r/ 2 ), and expressing the velocity v 
via the particle Lorentz factor 7 , Equation © can be 
presented in the following form 


sma 


me 2 7 


V 1 - l h 2 - 


(7) 


cos 2 a eB ro 

2 Note that at the motion with a constant absolute speed, the 
curvature radius by definition is equal to R c = v 2 /\a\. 


Apparently, for any given 7 this equation defines the an¬ 
gle a. Equation m indeed has solutions in the form of 
Equation ©» which describes the helicoidal motion de¬ 
pending on the parameters 7 and ro. Note that a motion 
with strictly circular trajectories is not possible. Thus 
the particle trajectory can be interpreted as a motion 
along the field line with a simultaneous drift in the per¬ 
pendicular direction. The positive and negative charges 
drift in opposite directions. 

The condition for a small step in the helix, \a\ -C 1, 
requires 


1« 7 « LAe = 6 x 10 « 


me* 


( B ) 

( r ° ) 

V10 11 G ) 

V 10 kmV 


( 8 ) 


Note that for typical pulsar parameters this condition is 
readily satisfied, and thus 


me 2 7 
eBr 0 ’ 


\a\ < 1. 


(9) 


Under the condition of Equation (j 8 ]), the curvature ra¬ 
dius R c of the trajectory of the particle located at the 
distance ro does not depend on the particle energy; it 
practically coincides with the curvature radius of the field 
line ro. Note that this conclusion is correct for any config¬ 
uration of the magnetic field. Indeed, the short segments 
of the field line can be treated as a circle, thus when satis¬ 
fying the condition of Equation ((HJ) the particle can move 
at this segment strictly along the drift trajectory. 

For the calculation of the prompt emission of the par¬ 
ticle moving along drift trajectory, the curvature of the 
field line can be used instead of the curvature of the tra¬ 
jectory because R c « ro- However, one cannot use the 
magnetic field line instead of the drift trajectory as a tra¬ 
jectory for the particle. As shown below, the characters 
of the motion of the particle initially moving along the 
magnetic field line and the drift trajectory are consider¬ 
ably different. Therefore, the interpretation of the cur¬ 
vature radiation as the radiation of the particle moving 
along the magnetic line is not correct. The strict formu¬ 
lation of the curvature radiation corresponds to the ra¬ 
diation of the particle moving along the drift trajectory. 
This is the definition that we will use for the curvature 
radiation throughout this paper. 

The intensity of the curvature radiation is defined as 


T 2 e 2 c 7 4 2 e 2 c 7 ' 

cu" = — R2~—- 


( 10 ) 


0 


When the particle moves close to, but not exactly along 
the drift trajectory, the intensity of radiation can devi¬ 
ate from Equation m- Therefore one should carefully 
investigate the trajectories in the vicinity of the drift tra¬ 
jectories. For this purpose, we invoke the Hamiltonian 
formalism. The vector potential can be taken in the form 
A x = 0, A y = 0, A z = A(r), therefore the azimuthal 
component of the field B{r) = — dA/dr. In the cylin¬ 
drical coordinates, the Hamilton function can be written 
as 


H = C yP r 2 + ^ + (P z ^^A(r)) 2 +m 2 c 2 , (11) 

where P r , and P z are the generalized momenta cor¬ 
responding to the coordinates r, <f> and z, respectively. 
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Because of the azimuthal symmetry of the magnetic field 
and its homogeneity along the z axis, <f> and z are cyclic 
coordinates. Therefore P ^ and P Zl as well as the en¬ 
ergy, are integrals of motion. The presence of three in¬ 
tegrals of motion allows us to reduce the problem to the 
one-dimensional case. Below we will limit the treatment 
to the case of ultrarelativistic particles and small pitch- 
angles. 

For given P ^ and P z , the radial motion can be consid¬ 
ered, as it follows from Equation m, as a motion in the 
field with the effective potential: 


U eS (r) = ^- + (P z ~-A(r)) 2 . (12) 

c 

The first and second derivatives of this potential are 


2P^ o 

K s (r) =- 1 + —B(r)(P z - - A(r)) , (13) 

r 3 c c 

and 

UM = & Jf + ^B\r){P z - e -A{r)) + %B\r). (14) 
r* c c c z 


Assuming that at r 
can write 


P* - 


= Tq, the derivative f/' ff (r 0 ) = 0 , we 


-AM 

c 


C P 2 

r%eB(r 0 ) ’ 


(15) 


t/'ffM 


2P | f 0 , r 0 B'(r 0 ) \ 

r o V B(r 0 ) J 


+ ^B 2 (r 0 ). (16) 

c 2 


If B decreases with r as a power-law, B ~ 1 /r* 5 , then 
for S < 3 both terms in Equation (fTTH) are positive, and 
therefore U e g has a minimum at ro- Furthermore, by 
assuming that the condition of Equation © is fulfilled, 
we can keep only the last term in Equation (1161) . 

For integrating the equations of motion, we apply the 
standard method that is used in the theory of small oscil¬ 
lations. In the vicinity of the minimum point, the func¬ 
tion U e g(r) can be approximated by a parabola: 


Ueg(r) = Pq + m 2 w 2 (r - r 0 ) 2 , 


where 


1 

2 m 2 


U? s (r 0 ) 


eB(rp) \ 2 
me ) 


(17) 

(18) 


Po = t/eff (r 0 ) 




n y 

rg?n 2 w^ J 


(19) 


Here oj c is the cyclotron frequency. 

In these denotations, the Hamilton function can be 
written in the form 

H = cyj P 2 + m 2 uj 2 (r — r 0 ) 2 + P 0 2 + m 2 c 2 , (20) 

and correspondingly the equations of motion are 


where E is the particle energy. The solution of these 
Equations gives 

D . W c t /OO'l 

r = ro — p cos-, P r = mpuj c sin-, (22) 

7 7 

where 7 = E/(mc 2 ) is the Lorentz factor of the particle, 
and p is an arbitrary const ant 0 related to the particle 
energy through the relation 


E = 



+ m 2 p 2 uj 2 + m 2 c 2 . 


(23) 


The condition of the applicability of Equation (1221) is the 
smallness of p compared to ro, as well as to the character¬ 
istic linear scale on which the magnetic field is changed 
significantly. 

Now we should find the time-dependencies of other co¬ 
ordinates. We note that in this region 


= 2 (r - r 0 ) 

dP,f, my r 2 7717 7 q \ ro 

Using Equation (l22l) . we find 


Substituting in 


the following approximate expression for A(r), 


■ (24) 


AI, 

7 

2 p uj c t\ 

— cos- , 

(25) 

777,7 r o 

V 

^0 

7 J 


7 

2^7 

. W c t\ 
sin- . 

(26) 

777,7 r o 

V 

Uc'1‘0 

7 j 

II 

cal qj 

11 

c 2 
' E 

( P ’~ 

- ^A(r)) 

(27) 


A(r) = A(r 0 + (r - r 0 )) « A(r 0 ) - B(r 0 ) • (r - r 0 ), (28) 
and using Equations (fl5l) and (l22l) . we find 


i = 


2 = 



7 w c m 2 rg 
P 2 


7 w c m 2 rg 


puj c u c t 
-cos-, 


7 7 


LU c t 

t — p sin-. 

7 


(29) 

(30) 


In order to make the analytical expressions more con¬ 
venient to work with it is useful to use, instead of P^ and 
p, the following parameters 




me'yro ’ 


p± 


PU C 
7 c 


(31) 


and introduce one more parameter, 


Pd 



C7W c 777 2 rg 


CT^f _ 

roU) c ‘ 


(32) 


Since P^ is the projection of the particle momentum 
on the axis z, /3|| is the parallel component of the velocity 
relative to the magnetic field, while (3 ± is the component 
of the velocity perpendicular to the drift trajectory, and 
is the drift velocity (all in units of c). 


_ dH _ c 2 P r ■ _ dH _ 
~~dP~ r ~ “E - ’ r “ ~ 


2 2 2 
E 


(r - r 0 ), 
( 21 ) 


3 Note that here the origin of the coordinate system is arbitrary, 
therefore instead of t one can use t — to, where to is the second 
arbitrary constant. Below we omit the arbitrary constants that 
can be added to </> and z. 
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For the new variables we find 

r/c = /3 ± sinr , i/c = /3 D — /3 ± cost , 


Ai c 


r o 


1 + 


2 Pj Pd 


COST 


(33) 

(34) 


where r = w c t/7, and the azimuthal component of the 
velocity is 


= r<j> = c/3|| 1 + 


Al Pd 


COST 


(35) 


For (3j_ = 0, this solution coincides with the previously 
obtained exact solution given by Equation ©, which de¬ 
scribes the motion along a drift trajectory (in this case, 
along the helicoidal line). Note that in this case the quan¬ 
tity /3 d plays the role of v z and describes the relativistic 
curvature drift. The same trajectory can be found if we 
average Equations ( 1351 ) - ( 1551 ) over time. Then this solu¬ 
tion can be interpreted as the motion along the helix 
around the drift trajectory. In the general case of an 
arbitrary magnetic field, motion strictly along drift tra¬ 
jectory is possible only locally, as long as the curvature of 
the field line is approximately constant. At /3 ± 7^ 0, the 
solution is approximate; it is correct when the following 
condition is fulfilled: 


P± ‘C /?||, /3 D <C /3y . (36) 

The relation between /3j_ and (3 D can be arbitrary. Since 
/3|| « 1, the second condition in Equation (l36ll is equiva¬ 
lent to Equation ([5]). The same requirement can be for¬ 
mulated in terms of the smallness of the Larmor radius 
compared to the curvature radius of the field line. 

Note that Equation (l28l) does not include the term cor¬ 
responding to the so-called gradient drift. Formally, the 
gradient drift could be taken into account if we add to 
Equation (1281) the next (quadratic) term of expansion. 
However our study is essentially linear, so the inclusion 
of the gradient drift could hardly be justified because it 
would exceed the accuracy of the approach. 

For the chosen new variables, the velocity and its 
square are 

v = re r + + ze z , (37) 

and 

v 2 = c 2 (pl + Pl + pl) . (38) 

It is important to note that v 2 does not depend on t. 
In order to find the acceleration, we use the following 
equation 

a = 6 (v x B) = — (re- — ze r ), (39) 

mc'y 7 

where it is taken into account that the field has only one 
(azimuthal) component. The square of acceleration is 

a 2 = Oq (l — 2rj cos r + rj 2 ) , (40) 

where 

a 0 = c 2 /3^/ro and 7/ = /3±/p D . (41) 


This equation can also be derived through the consider¬ 
ation of kinematics. Indeed, taking into account in the 
equation 

a = r = (r - r</> 2 )e r + (2 r<j) + r^)e^ + ze z . (42) 

the smallness of the ratios Pj_/(3^ and P D /f3 ,1, we arrive 
at an expression identical to Equation (l40l) . 

For /3j_ = 0 we have a 2 = a 2 = const. This is the 
case when the particle moves strictly along the drift 
trajectory. For (3 ± = f3 D , the square of acceleration 
a 2 = 2oq(1 — cost). For time instants, when the ac¬ 
celeration becomes zero, the particle velocity is parallel 
to the field line. In the limit of (3 ± (3 D , the accelera¬ 

tion becomes independent of the time and the radius of 
field curvature: 


2 2 Pi 

a =a 0 w 

Pd 


c/3j u] t 


(43) 


It is well known (see e.g. JLandau fe Lifshitzl (119751 11 that 
this equation also describes the particle acceleration in 
the homogeneous magnetic field. While for applicability 
of Equation ( 1551 ) the fulfillment of the condition /3j_ 

/3d is sufficient, the relation between /3 ± and (3^ can be 
arbitrary. 

For an ultrarelativistic particle, the acceleration and 
the curvature radius of the trajectory are linked with 
a simple relation |a| = c 2 /R c , which at v = c, is the 
definition of the curvature radius. Thus, the curvature 
radius of the trajectory is 


R c {t) = r 0 (l - 277COST + rj 2 ) 1/2 , (44) 

where /3|| is set to be equal to 1, given that 1 — /3|| <C 1. 
The pitch angle a of the particle also depends on time: 

a{t) = (3 d (l - 277 cos t + rj 2 ) 1 ^ 2 , (45) 

and varies in the range 

I Pd Pi. I < a < Pd + Pd ■ (46) 

Note that the angle between drift trajectory and particle 
velocity is constant and equal to /3±. 

3. INTENSITY AND ENERGY SPECTRUM OF RADIATION 

In this section we compute, within the framework of 
classical electrodynamics, the intensity and the energy 
spectrum of the radiation of a charged ultrarelativistic 
particle with the perpendicular component of velocity 
P± ~ Pd when both the curvature of the drift trajec¬ 
tory and the curvature of the gyration are important. If 
the radius of the curvature of the trajectory R c ~ ?~o, the 
characteristic energy of the emitted photon is ~ hc~/ 3 /ro. 
The requirement of smallness of this energy compared to 
the energy of the particle mc 2r ) (the condition of applica¬ 
bility of classical electrodynamics) implies S7 2 <C mcrg. 
We assume that the perpendicular (relative to the drift 
trajectory) momentum is p ± <C py, but, at the same 
time, p ± me. Using Equation (l32l) . the latter inequal¬ 
ity can be presented in the form C7 2 /g> w c ro. This gives 
the following constraint on the Lorentz factor: 


w c r 0 /c <C 7 2 « mcr 0 / ft. 


(47) 
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The condition of the smallness of the lower limit com¬ 
pared to the upper limit, implies tko c <C me 2 (i.e., 
the magnetic field should not exceed the critical field, 
B cr = 2m 2 c 3 /3eh « 2.94 • 10 13 G). 

The energy radiate d away by the particle during a 
unite time interval is (lLandau fc Lifshitzlll975 ') 


I = 


dE 

dt 


2e 2 7 6 

3c 3 


(v x a) 2 /c 2 ) 


2e 2 c y 4 

— R 2 ' 


(48) 


Substituting R c from Equation (l44ll . we obtain 


I = 


2e 2 c 7 4 


(l — 2?7 COSt + r] 2 ) . 


(49) 


The intensity averaged over time is 

U > = ^(1 + ^ 2 )- ( 50 ) 

D f Q 


This formula can be written in as (I) = J curv (l + ?? 2 ), 
where / CU rv is defined in Equation ([Toll . Thus, for a given 
Lorentz factor the particle produces radiation with min¬ 
imum intensity during the curvature radiation regime. 

In the considered scenario, we deal with three charac¬ 
teristic times: the time of energy losses (cooling time), 
tcooi ~ mc7'o/(e 2 7 3 ), the period of oscillations, T ~ 7/w c , 
and the characteristic time of formation of radiation, 
At ~ r 0 /(c7). It is easy to show, using Equation (H71) . 
that At <C tcooi- Note that this condition is always sat¬ 
isfied in the framework of classical electrodynamics. The 
ratio A t/T ~ w c ro/(c7 2 ) also is small compared to 1, 
as it follows from Equation (1471) . Finally, for the ratio 
r/tcooi ~ e 2 7 4 /(TncWcrg), using the lower and upper lim¬ 
its of the Lorentz factor from Equation (1471) . we obtain 


Hujr 

*f -2 

1 mc z 


< 


T me 

< a 


^cool 


t— - (51) 

lUjJn 


where a^ = e 2 /he ~ 1/137 is the fine-structure constant. 

The smallness of At significantly simplifies the calcu¬ 
lations of the radiation spectrum. This allows us to ig¬ 
nore the changes of the particle energy and the curva¬ 
ture radius with time, and perform calculations at fixed 
prompt values for these parameters: E{t) = mc 2r y(t) and 
R c {t). Then the spectral flux density of the magnetic 
bremsstrahlung integrated over all the emission angles 
is described by the well-known expr ession for th e syn ¬ 
chrotron reg i me of radiation (see, e.g. iSchwingeil (Il949fl : 
IB aver et all (11973IB : 


Here 


and 


P(w,t) 


V3e 2 7 / uj_\ 

2n R c \f2*/ 



pOO 

F(x)=x K 5/3 (x')dx', 

J X 


(52) 

(53) 

(54) 


where K 5 / 3 is a modified Bessel function. The par¬ 
ticle energy loss rate is determined to be dE/dt = 
— f™P(u,t) doj; the calculation of this integral results, 
as expected, in Equation (1481) . 


Remarkably, in all cases under consideration, the shape 
of the energy spectrum is defined by the same function 
F(x). The relevant parameters only change the position 
of the maximum and the intensity. However, the function 
F(x) shou ld b e replaced by its quantum analogue (see 
Equation m) if the parameter \ = Besma/B cr > 1, 
where e is the photon energy in units of m e c , and a 
is the angle between the photon and the magnetic field 
(jBeres tetskii ct al. 1982]). Note that at such conditions 
the energy of the produced photon is close to the en¬ 
ergy of the radiating electron. The electron-positron pair 
production by a gamma-ray photon in the strong mag¬ 
netic field occurs when % > 1. In the curved magnetic 
field the angle between the photon and the magnetic field 
could become sufficiently large for production of electron- 
positron pairs. This could lead to the development of an 
electromagnetic cascade, provided that the optical depth 
is large. The cascade results in the formation of radia¬ 
tion, the spectrum of which is considerably different from 
that of the initial (synchrotron) radiation. 

The limits of the applicability of Equation (15^1) , which 
describes the energy spectra of the synchrotron and cur¬ 
vature radiation compone nts, are det ermined by the ap¬ 
proach proposed bv ISchwingeil (Il949f) . The radiation of 
the ultrarelativistic particle is concentrated in a narrow 
cone with the opening angle ~ I/7, and therefore is col¬ 
lected when the angle between the velocity and the di¬ 
rection of the observation is of t he o rder o f the same as 
~ I/7. The method of iSchwinged (11949H is based on 
the expansion of the trajectory in a small time inter¬ 
val. During this time interval the entire observable radi¬ 
ation should be collected. This means that the approach 
works if the particle velocity changes direction at an an¬ 
gle larger than I/7 during the time that the expansion 
of the trajectory is valid. The analysis of the local tra¬ 
jectory given by Equations (l33l)(l35l) gives the following 
limits of applicability in the curved magnetic field: 


P± > -- 



7 

1 

» 


7 


if P.L 

if Pi .<Pd- 


(55) 


The first condition corresponds to the case of syn¬ 
chrotron radiation and states that the perpendicular 
motion should be relativistic as is expected from the 
consideration of radiation in the homogeneous magnetic 
field. The second condition corresponds to the situation 
of the synchro-curvature radiation of a particle with a 
small perpendicular momentum. Small gyrations around 
the drift trajectory do not influence the applicability of 
Equation (1521) . In the limit (3± = 0, this condition simply 
implies that the motion should be relativistic. 

Thus if the requirements of Equations (1551) are fulfilled 
for the synchrotron and curvature radiation regimes, the 
shape of the prompt energy spectrum remains the same. 
The distinct feature of the synchro-curvature regime is 
the significant time variability of the curvature R c . We 
assume that the condition T <C t C ooi is satisfied (note 
that this condition does not contradict Equation (l5lT) j. 
Then, since the change of 7 during the period T is small, 
we can introduce the parameters averaged over time. 
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This was already done at the changeover from Equation 
dH) to Equation assuming that during period T 

the Lorentz factor 7 = const. 

In accordance with Equation (l44ll . 

Rc = r 0 /q(v, r ) , q{r],T) = {l-2r]COST+r] 2 ) 1/2 . (56) 


Here a new parameter q determines the difference be¬ 
tween the curvature of the particle trajectory and the 
curvature of the magnetic field line, and thus charac¬ 
terizes the deviation of the particle radiation from the 
curvature radiation. 

Now the averaged spectral flux density can be pre¬ 
sented in the form 


A 7T Tq y id-. 


where 


UJ 


q{rj,r)F 

\U*J Wo 


UJ 


w*g(?b r ) 


(57) 


dr , (58) 


and w* = 3 c 7 3 /( 2 ro). G is a function of two variables: 
x = w/w* and 77 = /3 ± / )3 D . If in Equation (l58l) we repre¬ 
sent the function F(x) in the form of Equation (1M1) . and 
change the order of integration, the integral over dr can 
be calculated analytically. This allows us to express the 
function G(x) in the form of a single integral: 


n OO 

G(x) = x / K 5 /s(x') ^{x/x') dx ', (59) 

Jx/{\+1)) 


where 

T 


(?)-{r 


arccos 


(x/x') 2 — 1—rj 2 

2 V 


§r<\l-V\- 


(60) 


The function G for different values of is shown in Fig¬ 
ure |H At 77 = 0 (the dotdashed line) we have a nominal 
curvature radiation. The curves correspond to the nu¬ 
merical integration of Equation (1581) , using for the func¬ 
tion F(x ) the analytical approximation given by Equa¬ 
tion (IBS ll . In two limits, x 1 and x <C 1, the function 
F(x) can be expressed as 


F(x) w 2 2 / 3 Tx 1/3 , i«l. (61) 

F (. x )~ \j^Y e ~ x , x>l. (62) 

Therefore at x <C 1, the function G(x ) differs from F(x) 
only by a factor that does not depend on x: 

G(x) = — [ (1 — 277 cost + 77 2 ) 1 / 3 dr x F(x ). (63) 

n Jo 

In the opposite limit, x 1, the integral can be calcu¬ 
lated using the standard saddle point method. This gives 
the following asymptotics: 

GM ^ (J 0- ap (~Tri)' <64) 

Formally, in this case the saddle point method can be 
applied if x (1 + r/) 3 / 77 . However, the numerical cal- 


gW&O 



Figure 1. Spectral flux density of the synchro-curvature ra¬ 
diation averaged over time. The solid lines correspond to 
numerical calculations performed for the values of the ra¬ 
tio 77 = /3±/0d = = 0, 0.5, 1, 1.5 and 2 (from top to bot¬ 
tom, respectively). The dotted lines are calculated using the 
asymptotic presentation of the function G(x) given by Equa¬ 
tion (1611) . The dot-dashed line corresponds to the curvature 
radiation (77 = 0 ) when the particle moves along the drift 
trajectory. 


dilations show that in the interval 0.2 < 77 < 3 the sim¬ 
ple analytical function given by Equation (1641) already 
provides an accuracy better than 10 % at x « 2 (see 
Figure [T|). 

The spectrum described by Equations <53 and (1581) 
for a given magnetic field with curvature radius ro de¬ 
pends on the value 77 and Lorentz factor 7 . In a course 
of motion and radiation these parameters change, modi¬ 
fying the spectrum. In numerical calculations described 
further it is more convenient to use the q-parameter and 
Lorentz factors for the determination of the spectrum. 
As it follows from Equation (1581) the peak of the dis¬ 
tribution shifts with the change of parameters as gy 3 . 
This means that the peak of the spectra at a given en¬ 
ergy depends on the regime of the radiation determined 
by q. 

The local form of q given by Equation (1561) is useful 
for the theoretical studies. In particular it is seen that q 
oscillates between |1 — 77 I < <7 < I +77 with a frequency 
7 times smaller than the cyclotron frequency. With loss 
of energy the oscillation becomes faster. For 77 < 1, the 
g-parameter oscillates around 1 , sweeping the bandwidth 
277 . At 77 > 1 the oscillations occurs around 77 with an 
amplitude of 2. Taking into account Equation (l40l) the 
g-parameter can be written as 


q = 


a 

a 0 


evB sin a 
mcj 



sin a 


(65) 


where a is the pitch angle, v is the particle velocity, and 
v 11 is the velocity along the magnetic field. The approxi¬ 
mation sign shows that we have used the condition v ~ c. 
Moreover, this expression assumes that the conditions 
given by Equation (1361) are fulfilled. This form is par¬ 
ticularly convenient in numerical calculations possessing 














all the properties of the local form, as will be confirmed 
below. The value of the g-parameter allows us to follow 
the evolution of the radiation regime. If radiation pro¬ 
ceeds in the curvature regime, the g-parameter is close 
to 1 with small variations (77 -C 1) around it. When 
the amplitude of oscillations is of the order of the av¬ 
erage value of q, the radiation proceeds in the synchro¬ 
curvature regime. When the average value of q is much 
larger than the amplitude (Ag = 2) of oscillations, the 
particle radiates in the synchrotron regime, where the 
g-parameter could be very large. Correspondingly, the 
radiation spectra are shifted to higher energies. An im¬ 
portant conclusion that follows from the consideration 
of g-parameter is that the spectra produced in the cur¬ 
vature radiation regime (g « 1 ) is shifted to the lower 
energies in comparison with the spectra produced in the 
synchrotron and synchro-curvature regimes. 

For the derivation of the average spectral distribution 
we can use also a different approach. Equation ( 1551 ) can 
be presented in the form (a -C 1 ) 


^/?d\ 

) 

( 66 ) 

After trivial integration over da, we obtain Equa¬ 
tion ( 1551 ) . The function 

p(a) = -[ drS(a-P D q{p,r)) (67) 

7T JO 

describes the distribution of particles over pitch angles, 
thus p{a) da is the probability of the pitch angle being 
between (a, a + da). Therefore Equation (1661) can be 
written in the form 


a {^) = U dT l 


G 



/ da—p{a)F 
Jo Pd 


uj*a ) 


( 68 ) 


As a result G can be interpreted as the function F aver¬ 
aged over distribution of pitch angles. Note that in the 
derivation of Equation (1681) the form of the function q is 
not specified, so Equations (l 68 l) and (l58l) are equivalent. 

The determination of the spectral distribution of cur¬ 
vature radiation through integ ratio n over the pitch an¬ 
gles has been used bv lEnsteinl (1197311 (see Equation (20) 
in that paper). However, one has to be cautious when 
applying this method. The pitch angle of particles mov¬ 
ing in the field with curved lines depends on time (con¬ 
trary to the case of motion in the homogeneous magnetic 
field). Therefore the distribution p{a) can not be a dopted 
in an arbitrary manner, as was done bv lEnsteinl ( 19731 ). 
This function should be derived from Equation (1671) or 
obtained from the solution of relevant kinetic equations. 
For example if q{p, t ) is given by Equation (l56l) . we ob¬ 
tain the following expression: 


p(a) = - 2a = (69) 

(Pbiv + ~ « 2 ) (« 2 -/%(»?- x ) 2 ) 

in the range of pitch angles fi D \p — 1 | < a < P D {p + 1 ), 
and p(a) = 0 outside that interval. Note that at rj ^ 1, 
we do not have particles with a null pitch angle. The 
integration of Equation (1681) over this distribution func¬ 
tion gives the same result as that obtained by integration 


over the time. 

Let us consider now the case when the accelerated elec¬ 
trons are not strictly unidirectional but are distributed 
within a narrow cone with an opening angle At/> <C 1 
and axis along the drift trajectory. We assume that all 
electrons have the same energy. For the standard syn¬ 
chrotron radiation, a slight change of the pitch angle does 
not affect the spectrum of radiation. However, in the case 
of the synchro-curvature radiation even a tiny spread of 
the angles relative to the drift trajectory may result in 
a significant change of the spectrum. It is convenient to 
express the angular distribution of particles through the 
parameter 77 , which is linked to the the angle between the 
velocity and the drift trajectory with a simple relation: 

= 77 / 3 d // 3 ||. Note that at p < 1, the angle 1. For 
calculations, the angular distribution of particles should 
be specified. Below we consider a Gaussian type distri¬ 
bution: 


g( 77 ) dp ~ e v dp , (70) 

with the mean value of the square of the perpendicular 
velocity 

</3i)=C Pi- (71) 

For the determination of the energy spectrum we have to 
average Equation (1581) over the time r and the parameter 
77 , i.e., calculate the integral 


(GW, = 1 (72) 

The energy spectra calculated for several values of £ 
are shown in Figure [2] The case of £ = 0 corresponds to 
the spectrum of the curvature radiation. 

As above, using the saddle point method, we can find 
the asymptotics for ( G) at large x. Rather simple al¬ 
though cumbersome computations, which we omit here, 
lead to the following result 

< G(x) ) ~ \f^73 C 1 "' 3 + C“ 1/S ) 2 e “" ' P3) 

where 

u= 3 c -VV /3 _ C - 2/ V /3 + _L. (74 ) 

This equation is derived under the following conditions: 
i>( and x 1 /£. However, the comparison of Equa¬ 
tion (1741) with accurate numerical calculations shows (see 
Figure[ll) that for ( ~ 1 this analytical presentation gives 
correct results at x > 0.5. 

Although the assumed Gaussian type angular distri¬ 
bution of electrons seems to be a quite reasonable and 
natural choice in the considered scenario, it would be 
interesting to investigate the dependence of the radia¬ 
tion spectrum on the specific angular distribution of the 
electron beam. In particular, in Figure (3J we show the 
radiation spectra calculated for the uniform distribution 
of electrons within the fixed opening angles of the beam. 
The solid curves in Figure [3] are obtained for the same 
values of the mean square of the perpendicular velocity 
as in Figure [2] In Figure [3] we show the asymptotic so¬ 
lutions given by Equation (l73l) for the spectra calculated 
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Figure 2. Spectral flux density of the synchro-curvature ra¬ 
diation averaged over the time, and the angle between the 
velocity and the drift trajectory, assuming a Gaussian type 
distribution given by Equation Cm]). The solid lines corre¬ 
spond to different values of £ = {/3±)//3 1, 3, and 5 (from 
bottom to top, respectively). The dotted lines are calculated 
using the asymptotic analytical presentation given by Equa¬ 
tion ([73]). The dashed-dotted line corresponds to the curva¬ 
ture radiation (£ = 0). 


Table 1 

The spectral index s at different points x = lj/u* of the spectrum 
( G(x)) given by Eq. (fT3l) for various values of £ = {/!) ) /d/j - 



2.0 

3.0 

5.0 

10. 

20. 

30. 

0 

1.65 

2.62 

4.59 

9.56 

19.5 

29.5 

0.5 

1.17 

1.67 

2.56 

4.45 

7.56 

10.2 

1.0 

0.96 

1.37 

2.10 

3.63 

6.14 

8.28 

2.0 

0.75 

1.08 

1.67 

2.91 

4.93 

6.65 

3.0 

0.64 

0.93 

1.45 

2.54 

4.31 

5.82 


for the Gaussian type angular distribution (the same dot¬ 
ted curves from Figure [2]). Figure [3] demonstrates that 
the value of the mean square of the perpendicular com¬ 
ponent of velocity describes the energy spectrum of ra¬ 
diation quite well, independent of details of the specific 
small pitch-angle distribution of electrons. 

Note that the spectra of the synchro-curvature radia¬ 
tion do not contain the characteristic exponential term 
e~°°, as is the case of the synchrotron or curvature ra¬ 
diation, but show a harder behavior. This can be seen 
from the asymptotics in Equation (l74l) . To demonstrate 
the tendency of the steepening of the energy spectrum 
beyond the maximum at x ~ 1, in Table Q] we show the 
evolution of the local slope of the spectral flux density 
(the spectral index s), which is determined to be 

fl = -d(log<G))/d(logw). (75) 

It is seen that whereas in the case of curvature radiation 
(£ = 0) the spectrum becomes very steep at x > 3, in the 
synchro-curvature radiation regime with £ = 2 or £ = 3 
a relatively hard spectrum can extend up to x ~ 10. 

4. NUMERICAL IMPLEMENTATION 


<gW&0> 



Figure 3. Spectral flux density of the synchro-curvature ra¬ 
diation averaged over the time and angle between the velocity 
and the drift trajectory, assuming a uniform angular distri¬ 
bution of electrons within the opening angle of the beam. 
The solid lines correspond to different values of the ratio 
£ = {0±)/P d'- lj 3, and 5 (from bottom to top, respectively). 
The dotted lines are the same as in Figure [2] They are calcu¬ 
lated for the Gaussian type of angular distribution using the 
asymptotic analytical presentation given by Equation <113. 
The dashed-dotted line corresponds to the curvature radia¬ 
tion (<( = 0). 


The analytical approach described above allows us to 
study the local properties of the particle trajectory and 
its radiation in the curved magnetic field. To solve the 
problem in the general case, namely, to take into account 
the energy losses of the particle during its propagation, 
we performed numerical integration of the equations of 
motion. The radiation properties have been studied for 
the dipole magnetic field. This is a good approximation 
for the case of a strong magnetic field in the vicinit y of 
compact astrophysical objects. As it is shown in lDeutschl 
(11955D . the instantaneous configuration of the magnetic 
field in the vicinity of a rotating star appears as a sta¬ 
tionary magnetic dipole. Furthermore, in Secti on [6] we 
use the approximation given by iDeutschl ( 1955 ) for the 
electric field in the vicinity of the star. The dipole mag¬ 
netic field has two distinct features. The first is the fast 
decrease of its strength with the distance, B(r) oc 1/r 3 , 
which has a strong impact on the radiation intensity. The 
second is the significant variation of the curvature of the 
field lines, with a change of the polar angle 6 from the 
dipole axis oc sin0/r. Therefore the radiation spectra in 
the vicinity of the pole and the equator have different 
characters. 

For the calculation of the exact trajectory of a particle 
in the curved magnetic field, we have used the equations 
of motion in the ultrarelativistic limit. In this limit the 
velocity has the constant value of the speed of light and 
changes only its direction. The radiation reaction force 
is directed opposite to the velocity. Thus the equation of 
motion takes the form 

me {^~dt +7 ^) = X ~ ^ 
























10 


where (3 is the velocity in units of c with \(3\ ss 1, B, 
B is t he magnetic field, and f is the radiation reaction 
force ijLandau fc Lifshitzll 19751 1. Here, the left side of the 
equation is expanded for the convenience of further trans¬ 
formations. Taking the scalar product of Equation (ITTH) 
with velocity f3 and using the fact that (3 2 = 1 and, 
therefore, /3^ = 0 , we obtain the following differential 
equation for the Lorentz factor 


= _]/[ 

dt me 


(77) 


Substitution of this equality back to Equation d76l) gives 
us 


d/3 

dt 


e 

mc'y 


(0 x B), 


(78) 


which has the same form as if we considered the motion 
without energy losses. However, we now have two equa¬ 
tions where the Lorentz factor is variable: Equations JZ3 
and (ZED- For the sake of convenience of the numerical 
treatment and comprehension of the structure of the sys¬ 
tem of equations, these equations (plus the equation for 
coordinates) can be written in the dimensionless form 
(see also Appendix lAl) : 


^r = «u9, 

(79) 

f? = £09 x b ), 

(80) 

7 1 = u 2 {(3 x b) 2 . 
dr 

(81) 

The system of equations depends on 

two dimensionless 


parameters 


Ml 


me 2 70 
eBpRp 


and 



e 3 B 0 2 
m 2 c 4 ^ 0 ’ 


(82) 


where 70 is the initial Lorentz factor, Rp is the charac¬ 
teristic distance to the radiating region from the dipole 
(the distance to the outer gap or polar cap), B 0 = 
B * (R*/Rp) 3 is the characteristic magnetic field at this 
distance with P* and P* being the star radius and the 
magnetic field at the magnetic pole of the star, m is a 
particle mass, and c is the speed of light. Here we have in¬ 
troduced the following dimensionless variables: r' is the 
coordinate in the units of Rp , (3 is the velocity of the par¬ 
ticle in the units of c, £ = 70 /”f is the ratio of the initial 
and current value of the Lorentz factor, r = teBp/mc'yp 
is the characteristic time in units of the initial gyration 
period, and b is the dimensionless dipole magnetic field 
that is expressed as 


(83) 

where /x is the unit vector in the direction of the dipole 
axis, and n = r’ /r’ is the unit vector in the direction of 
the particle position. 

It should be noted that, depending on the specific con¬ 
ditions characterizing an astrophysical source, the pa¬ 
rameters Mi and U 2 may differ from each other by many 
orders of magnitude. For instance, for typical parame¬ 
ters of the polar cap model of the pulsar magnetosphere 
Bp = 10 12 G, Rp = 10 6 cm, and 70 = 10 8 , we have 


Mi ss 1.7 x 1CH' and M 2 « 1.1 x 10 12 . Thus, we deal with 
a non-stiff problem. The implicit Rosenbrock method 
has been used for integration of this system of differen¬ 
tial equations. 

The calculations of the trajectory were performed for 
different initial conditions. The initial position is deter¬ 
mined by the radius Rp and the polar angle Op relative 
to the magnetic dipole axis. The radiation spectrum was 
calculated for different initial pitch angles ap and the 
initial Lorentz factors 70 . The information regarding the 
particle trajectory and its energy allows us to derive the 
^-parameter from Equation (l65l) . and then calculate the 
radiation spectrum at any moment of time using Equa¬ 
tion (|58l) . The cumulative spectrum (i.e., the spectrum 
integrated over the time), can be calculated for any given 
initial parameters. 

Finally, in the case of a very strong magnetic field 
and/or very large Lorentz factor, the particle may radi¬ 
ate in the quantum regime. More specifically, when the 
parameter \ = B^/ sina/B cr > 1, Equation (ED s houl d 
be replaced by its quantum analogue (see Equation (IB6D 1 

§■ =- x 6)2 - 5 , (84) 

dT (l + m 3 ((/3 x b)/£) 2 / 3 )“ 

where M 3 = 1.07 x 10“ 9 (i?o7o) 2 ^ 3 - 

5. ASTROPHYSICAL IMPLICATIONS 

In this section we explore some possible realizations 
of the synchrotron and curvature regimes, as well as 
transitions between these two regimes of radiation (the 
synchro-curvature regime) in the context of two specific 
astrophysical scenarios. Namely, we discuss the radiation 
of electrons and positrons in the pulsar magnetosphere 
for the outer gap and polar cap models. 

Before proceeding, we first make some esti¬ 
mations of the basic parameters in the dipole 
magnetic field. The radius of light cylinder is 
Ric ~ 1.57 • 10 8 P CT . cm where P cr = P/33 ms and 
P is the period of pulsar. The strength of mag¬ 
netic field is B k 1.28 ■ 10 5 B 12 P~ 3 R~ 3 G, where 
B 12 = P/10 12 G and R = r/Ri c . The drift velocity is 
(3d ~ 1-26 ■ lO _ 3 77 i? 2 Hf 2 1 P 2 r sin 0 , where 77 = 7 /IO 7 , 
and 0 is the polar angle from dipole axis. Here we use 
the formula /3d = 'jmc 2 /eBrp, where rp is the curvature 
of the magnetic field lines and /3m =1 (compare with 
Equation (1371) 1. Then according to Equation © 
q ~ 0.79 ■ 10 3 a 7 ^ 1 P _ 2 Pi 2 P “ 2 sin -1 0, where a is the 
pitch angle. The characteristic energy of the curvature 
radiation is e curv ~ 2.8 ■ 10 8 7 3 P _ 1 P “ 1 sinf? eV. Thus, 
an electron with q = a/po, which determines the 
regime of the radiation emits photons with characteristic 
frequency ep = qe curv with the intensity I = q 2 I cur v , 
where I C urv is the intensity of the curvature radiation. 


5.1. Outer Gap 

The high-energy gamma radiation from pulsars is be¬ 
lieved to origin ate from the oute r gap of the pulsar 
magnetospher e (jCheng et alJ 119861 : iTakata et al J 120041 : 
Hirotani 2008). Here we present the results of our cal¬ 
culations of the radiation of electrons (positrons) in the 



















11 



Figure 4. Cumulative (integrated along trajectory) radiation spectra of electrons calculated for the outer gap model in the pulsar 
magnetosphere. The curves are obtained for the initial Lorentz factors of electrons 7 = 10® and 10 7 (left panel), and 7 = 5 X 10 7 and 10 8 
(right panel). The spectra calculated in the approximation of the motion strictly along magnetic field lines are shown by solid lines. Other 
spectra are calculated using the exact trajectories of electrons with different initial directions relative to the magnetic field lines: along the 
drift trajectory (dotted lines) (not seen on the left panel because of the coincidence with solid lines), along the magnetic field line (dashed 
lines), and at pitch angle a = 1 0/3 d in the meridional plane opposite to the normal vector of magnetic field lines (dashed-dotted lines). 



Figure 5. Time evolution of the q-parameter and the electron 
Lorentz factor in the outer gap model (complementary to Figure[4]). 
The four panels correspond to the initial Lorentz factor of electrons 
7 = 10®, 10 7 ,5 X 10 7 ,10 8 and their initial direction along the drift 
trajectory. 

dipole magnetic field at the location of the outer gap. 
The position of the outer gap is placed above the last 
open held line. The inner boundary of the gap is at the 
null surface where the GoldreichJulian charge density is 
zero, and the outer boundary is taken on the surface of 
the light cylinder. We adopt the parameters of the Crab 
pulsar: the radius of the star i?* = 10 6 cm, the rotation 
period P = 33.5 ms, and the magnetic held at the pole 
B* = 10 12 G. We consider a non-aligned pulsar with an 
angle of 45° between the rotation axis and the magnetic 
dipole axis. The initial particle position is set at the 
distance r; n i t = 0.5 Ri c from the null surface along the 
last open held line, where Ri c is the radius of the light 
cylinder. In numerical calculations, the particle is fol¬ 
lowed up to the intersection with the light cylinder. For 
parameters typical for the Crab pulsar, these conditions 


Figure 6. Same as in Figure [5] but for the initial direction of 
electrons along the magnetic field line. 

correspond to 8 0 = 51.5°, Rq = 1.2-10 8 cm, B 0 = 5.6-10 5 
G, and Ri c = 1.6 • 10 s cm, where 8q is the polar angle 
relative to the magnetic dipole axis. 

The energy spectra and the radiation regimes were 
studied for different initial directions and Lorentz fac¬ 
tors of electrons. The resulting spectra are shown in 
Figure |U The curves in Figures [5] [7] demonstrate the 
time-evolution of the g-parameter (blue lines, right scale) 
and the Lorentz factor normalized to its initial value. 
These complementary plots allow us to simultaneously 
watch the energy loss rate and the regime of the radia¬ 
tion characterized by the value of the g-parameter. The 
saltatory behavior of some curves for the g-parameter is 
produced by small number of data points in the region 
when the energy loss rate is low. 

It follows from the above discussion that the particle 
with a given energy radiates less intensively and pro¬ 
duces less energetic photons when it moves along drift 
trajectory (i.e., in the curvature radiation regime). In- 
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Figure 7. Same as in Figure |5] but for the initial direction of 
electrons at the angle 10 0d- 


deed, locally such motion corresponds to the case when 
7] <C 1 or /3j_ <C Pd in Equations ([50]) . (l56l) . and ([58]) . It 
means that q oscillates around 1 with a small amplitude 
of r). Then the peak of the spectrum given by Equa¬ 
tion (1581) does not shift from the position determined in 
curvature radiation. To check this conclusion, we set the 
initial direction of the particle along the drift trajectory, 
thus fixing the initial velocity at the angle (3d from the 
magnetic field line toward the binormal vector. The cor¬ 
responding cumulative (integrated along the trajectory) 
energy spectra of radiation are shown in Figure 0] The 
complementary plots are shown in Figure [5] We can see 
that, indeed, for the initial Lorentz factors 70 = 10 6 and 
70 = 10' the curvature radiation is less energetic com¬ 
pared to other cases and the rates of energy losses are 
minimum as well (compare with corresponding curves in 
Figures |6] and [T| . 

However, the situation is different for the larger ini¬ 
tial Lorentz factors 70 = 5 • 10 7 and 70 = 10 8 shown 
on the right panel of Figure [4] It is seen that for the 
initial direction along the drift trajectory, the radiation 
spectra extend to higher energies than in the case of the 
initial direction along the magnetic field. This can be 
explained by the very intense energy losses occurring be¬ 
fore the particle has made the first gyration (the first 
oscillation of q-parameter). During this time interval, 
q ss 1 for the initial direction along the drift trajectory 
and q ~ 0 for the initial direction along the magnetic 
field. Equations (l57l) and (l58l) show that larger values 
of q give a more energetic radiation. Correspondingly, 
more energetic radiation is produced in the case of the 
initial direction along the drift trajectory. However, after 
many gyrations the energy losses in the case of the initial 
direction along the drift trajectory become, as expected, 
less intense compared to the case of the initial direction 
along the magnetic field line. 

In Figure [I] we compare the exact calculations of radi¬ 
ation in the curvature regime with the results obtained 
under the assumption that the electron moves strictly 
along the magnetic field (as is routinely assumed in many 
papers on curvature radiation). As expected, for the 
relatively modest initial Lorentz factors 70 = 10 6 and 


70 = 10 7 , on the left panel of Figure [2 the curves 
for these cases merge. For larger initial Lorentz factors 
70 = 5 • 10 7 and 70 = 10 s the differences are not negligi¬ 
ble, but are still small. 

We call the attention of the reader to the regimes of 
radiation demonstrated by the curves in Figure [5] For 
70 = 10 6 and 70 = 10 7 , the q-parameter equals unity, 
indicating that the radiation proceeds in the curvature 
regime. However, this equality is not exact, as is demon¬ 
strated for 70 = 10 7 where q-parameter slightly oscillates 
around unity. These small oscillations correspond to the 
fine gyration around the drift trajectory. As discussed 
above (see Equation (1551) 1 the presence of such fine per¬ 
pendicular motions does not influence on the applicabil¬ 
ity limits. 

For 70 = 5 • 10 7 and 70 = 10 8 , the q-parameter has a 
more complex behavior. The increase at the beginning 
is defined mostly by the fast energy losses. The decrease 
is determined by the combination of several factors, such 
as the reduction of the magnetic field strength and the 
change of its curvature. The increase of q-parameter indi¬ 
cates that the radiation occurs in the synchro-curvature 
or the synchrotron regimes when the radiation due to 
curvature of the magnetic field line is less important (for 
70 = 5 • 10 7 ) or simply negligible (for 70 = 10 s ). Accord¬ 
ing to Equation (l58l) . the energy of the radiation maxi¬ 
mum scales as ~ 7 3 q. The q-parameter reaches the max¬ 
imum when a considerable amount of energy has been 
lost. Therefore, in spite of large q, the peak of radiation 
shifts toward low energies and does not affect the cu¬ 
mulative spectrum. The interesting feature can be seen 
at first moments when the energy oscillates with the q- 
parameter and the minimums of q-parameter correspond 
to flatter parts of the Lorentz factor evolution curve. 

The radiation spectra of electrons launched along the 
magnetic field line are slightly more energetic, except for 
70 = 5 • 10 7 and 70 = 10 8 as discussed above. Initially, 
the radiation for 70 = 10 6 and 70 = 10' is in the synchro¬ 
curvature regime (see Figure 0, although for 70 = 10 6 
the most of the energy is lost in the curvature regime, 
which occurs fast. Therefore the spectrum in this case 
almost coincides with the spectrum for the initial direc¬ 
tion along the drift trajectory. For 70 = 10 7 most of 
the energy is lost in the synchro-curvature regime, thus 
the spectrum is shifted to higher energies by a factor of 

q« 2. 

For illustration of the effect related to initial pitch an¬ 
gles larger than (3d, we show the case of the initial di¬ 
rection deflected at the angle 10/3 d from the magnetic 
field line toward the direction opposite the normal vec¬ 
tor. The corresponding spectra indicated in Figure [4] by 
dashed-dotted lines are shifted toward higher energies. 
Although the q-parameter (Figure [7]) reaches large val¬ 
ues, most of the energy is lost at initial stages at q ~ 10. 
Thus the spectra are shifted by q « 10 compared to the 
curvature radiation spectra. 

We should note that the energy spectra of radiation 
produced in all regimes at highest energies contain an 
exponential cut-off similar to the spectrum of the small- 
angle synchrotron radiation. However, since the radia¬ 
tion spectrum is very sensitive to the pitch-angle, a pop¬ 
ulation of electrons with similar energies but different 
angles can result in a superposition spectrum with a less 
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abrupt cut-off. The condition for the realization of such a 
spectrum is that the distribution of electrons over the an¬ 
gles around the drift trajectory should be wider than [3d- 

5.2. Polar Cap 

In the polar cap model the electron radiates in the 
region located close to the surface of the neutron star, 
where the magnetic field is much stronger than in the 
outer gap model, approaching to B m 10 12 G. This re¬ 
sults in much faster damping of the perpendicular com¬ 
ponent of motion. The very small drift velocity [3d im¬ 
plies that the drift trajectory and the magnetic field line 
almost coincide, thus even a small deflection from the 
magnetic held line produces radiation that is quite differ¬ 
ent from the curvature radiation. However, the transition 
to the curvature radiation regime occurs very fast. In the 
curvature regime electrons radiate more energetic pho¬ 
tons than in the outer gap model because of the stronger 
curvature of the held lines. The curvature of the dipole 
held lines behaves like ~ sin 9/r, and with the decrease of 
the radius by two orders of magnitude and the decrease 
of the polar angle by an order of magnitude (do ~ 1 °), 
the curvature is larger by an order of magnitude com¬ 
pared to the case of outer gap model. Correspondingly, 
the maximum energy of the curvature radiation in the 
polar cap is an order of magnitude higher than in the 
outer gap. 

The energy spectra of radiation calculated for the po¬ 
lar cap model are shown in Figure [ 8 ] The complemen¬ 
tary plots for the evolution of the ^-parameter and the 
electron Lorentz factor are presented in Figures [Mill 
The initial position of the particle is Rq = 10 6 cm, and 
do = 3°. The spectra indicated by solid lines correspond 
to the case when the initial direction of the particle is 
along the magnetic field line. In this case the radiation 
is in the curvature regime. Initially, the radiation shortly 
proceeds in the synchro-curvature regime, but abruptly 
turns to the regime with fast oscillations around q ss 1 
caused by fine gyrations (Figure [9]). 

In the case of an initial pitch angle ao = I /70 the 
abrupt change of regimes leads to an interesting feature 
in the cumulative spectra (dashed lines). Because of 
small changes of the energy and fast changes of the q- 
parameter, an energy spectrum is formed consisting of 
two peaks. The peak at higher energies is produced in 
the synchrotron regime (g> 1 , see Figure [TUI) , while the 
lower energy peak is due to radiation in the curvature 
regime. The double-peak structure disappears for large 
initial pitch angles. For example, for the pitch angle 
IOO /70 the transition to the curvature regime is very fast 
and the electron enters into this regime with dramatically 
reduced Lorentz factor. Thus the peak of the curvature 
radiation is not only shifted to smaller energies, but also 
too weak to be seen in the cumulative spectrunfl. 

In very strong magnetic fields, namely when the pa¬ 
rameter x = sma/B cr > 1 , the radiation is pro¬ 

duced in the quantum regime. Let’s assume that the 
initial pitch angle is inversely proportional to the initial 
Lorentz factor, ao = «/ 7 o- This makes the condition of 
radiation in the quantum regime independent of 70 : 

4 Note that the pitch angle IOO/ 70 , which we treat as large, is 
still extremely small, ~ 2 ( 70 / 10 7 )~ 1 arcsec. 



Figure 8. Cumulative (integrated along trajectory) radiation 
spectra of electrons calculated for the polar cap model in the pul¬ 
sar magnetosphere. The curves are obtained for different initial 
Lorentz factors of electrons 7 = 10®, 10 7 ,10 s , and for different ini¬ 
tial directions relative to the magnetic field lines: along the mag¬ 
netic field line (solid lines), and for two pitch angles I /70 (dashed 
lines) and IOO /70 (dot-dashed lines). 


B 1 - 10 12 G 

x = _ 78 ii lao = ____ a «,3.4. 1 °-«. (85) 

Thus, at the initial pitch angle with a > 30, the elec¬ 
trons radiate in the quantum regime. Dashed-dotted 
lines in Figure [ 8 ] present radiation spectra for the ini¬ 
tial pitch-angles ao = IOO/70. Note that in the quantum 
regime almost the entire energy of the parent electron is 
transferred to the radiated photon. Therefore we should 
expect an abrupt cutoff in the radiation spectra. This 
effect is clearly seen in Figure [ 8 ] (dot-dashed curves cor¬ 
responding to the initial pitch angle ao = IOO/70). The 
gamma-rays produced in the quantum regime are suf¬ 
ficiently energetic to be absorbed in the magnetic field 
through the e + e _ pair production. This leads to the de¬ 
velopment of an electromagnetic cascade in the magnetic 
field. The spectrum of cascade gamma-rays that escape 
the pulsar magnetosphere will be quite different from the 
spectra shown in Figure [H 

6 . MOTION AND RADIATION IN ELECTRIC AND 
MAGNETIC FIELDS 

To study a more realistic picture of the radiation, 
one should consistently determine the initial pitch angle. 
This can be done if the acceleration of the particle is con¬ 
sidered. First we consider the motion of a charged par¬ 
ticle in the crossed homogeneous electric and magnetic 
fields. If the fields are not perpendicular to each other, 
one can always find a reference frame in which they are 
parallel. In this reference frame the particle completely 
loses the momentum perpendicular to the common di¬ 
rection of the electric and magnetic fields. At the same 
time, the electric held infinitely accelerates the particle. 
Eventually the particle accelerates along the common di¬ 
rection of the electric and magnetic Helds. In an arbitrary 
reference frame the particle is accelerated along a recti¬ 
linear trajectory. However, the direction of the motion 
is 
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Figure 9. Time evolution of the ^-parameter and the electron 
Lorentz factor in the outer gap model (complementary to Fig¬ 
ure The three panels correspond to the initial Lorentz factor 
of electrons 7 = 10 6 , 10 7 , 10 8 and their initial direction along the 
magnetic field line. 



Figure 10. Same as in Figure [9] but for the initial direction of 
electrons at the pitch angle 1 / 70 - 




^+^4 <-> 


where 


Q = 


£ 2 +B 2 \(£ 2 B 2 


+ ( £B )2, (87) 


and + and — signs correspond to positively and neg¬ 
atively charged particles, respectively. Then the pitch 
angle can be expressed as 


sin a = 



Figure 11. Same as in Figure [9] but for the initial direction of 
electrons at the pitch angle 100 / 7 o- 

where a = S’/B is the ratio of the electric and magnetic 
fields, and 0 is the angle between them. Assuming that 
the electric field is typically smaller than the magnetic 
Held, we obtain 

a ~ Pe = cr sin 9 -C 1 (89) 

which implies that the pitch angle equals the electric 
drift velocity (in units of the speed of light). Thus, in 
the crossed homogeneous electric and magnetic fields, 
the asymptotic motion of the particle after damping of 
gyration is a constant acceleration along the rectilinear 
trajectory. The direction of the trajectory constitutes 
the angle /3g with the direction of the magnetic field. 
Note that due to this acceleration the particle radiates. 
However, the energy losses do not depend on the Lorentz 
factor and are negligibly small compared to the energy 
losses during the particle motion along the curved tra¬ 
jectory. 

In the curved electric and magnetic fields the accelera¬ 
tion of the particle is restricted by the energy losses due 
to the curved trajectory. Once the characteristic times of 
the acceleration and radiation become equal, the particle 
moves along a drift trajectory determined by the elec¬ 
tric and magnetic fields. During radiation the particle 
changes its direction of motion. At the same time it is ac¬ 
celerated by the electric field in the direction determined 
by Equation (l86l) . As in the case of crossed homoge¬ 
neous electric and magnetic fields, the component of the 
momentum perpendicular to the drift trajectory, which 
is responsible for the gyration of the particle, is damped 
due to radiation. Thus, in the equilibrium between accel¬ 
eration and radiation, the particle moves along drift tra¬ 
jectory practically without gyration with the pitch angle 
/3 e- If the electric or magnetic field changes slowly along 
trajectory, it gradually adjusts itself to the equilibrium 
state of the drift trajectory. However, if fields vary fast, 
it can lead to a strong radiation until the particle gets a 
new equilibrium drift trajectory. Such conditions might 
appear when the particle escapes the acceleration gap 
and the field-aligned (along magnetic field lines) electric 
field is screened. 

Further we study the radiation of the particle inside 
and outside the outer gap. We check how the character 
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of the radiation changes depending on the variation of 
the electric field. Specifically, we consider the case of 
screening the fieldaligned electric field when it decreases 
exponentially from the border of the outer gap. The 
geometry of the model is presented in Figure [T2] As 
before, the angle between the rotation and dipole axes is 
45°. The strength of the magnetic field at the pole of the 
star is B * = 10 12 G. The period of rotation of the star 
is P = 33 msec. The configuration of the electric and 
magne tic fields is given by the approximation of lDeutschl 
( 1955 1 for the fields in the vacuum in the vicinity of a 
rotating star (equations (18) and (19) in iDeutschl (|1955f) ). 
The magnetic field in this approximation is just a dipole 
magnetic field. The outer gap is located between the null 
surface and the last open field line. We set the border 
of the gap at Xb = 0.7 Ri c from the rotation axis (see the 
black vertical line in Figure H2l) . where Ri c is the radius 
of the light cylinder. In our calculations we assume that 
the field-aligned electric field is screened starting from Xb 
according to the exponential decay: 

^ =4||e-(^) , (90) 

where <£q|| is the original component of the electric field 
parallel to the magnetic field (as it would change with¬ 
out screening). Note that the perpendicular component 
stays unmodified. Changing the width of the screen¬ 
ing W, we can explore different models for the variation 
of the electric field. The three cases that we investi¬ 
gate, W = 10 ~ 3 Ric, 10 ~ 2 Ri c , and 10 ~ 1 Ri c , cover a large 
range of values of this parameter. The calculation of 
the trajectory starts from the null surface at a distance 
x = 0.4i?; c from the rotation axis, with initial Lorentz 
factor 70 = 100 and velocity directed along magnetic 
field line. 

Together with calculations using the exact equations of 
motion given by Equation m, we performed the cal¬ 
culations in the drift approximation presented by Equa¬ 
tion (ESI). The approach of the jdrift appro ximation has 
been used in the work iKalapotharakos et al. (120121 ) to 
calculate the light curves for the different models of the 
magnetosphere. Here it is important for us to inves¬ 
tigate how well this approach reproduces the radiation 
spectra. Apart from the described configuration of the 
electric field, we carried out the calculations for more 
simple electric fields when their strength is proportional 
to the strength of the magnetic field and their direction 
is always at a constant angle to the direction of the mag¬ 
netic field at a given point. The results for a broad range 
of values for the angles and strength ratios have shown 
a very good agreement between the exact and approxi¬ 
mate approaches. We do not show them because the case 
of the configuration considered here is quite representa¬ 
tive. The explanation of the good agreement between 
radiation spectra calculated in exact and approximate 
approaches is that in both cases the curvatures of the 
trajectories are very close. As discussed previously, in 
the equilibrium between the acceleration and radiation 
the particle basically moves along drift trajectory, which 
is assumed a priori in the drift approximation. 

The projection of the trajectories on the XZ plane cal¬ 
culated in both approaches are presented in the Fig¬ 
ure [T2] The trajectories outside the gap are shown 


for the case of the slow variation of the electric field 
W = 10 ~ 1 Ri c when the scale of change is of the same 
order as the scale of the gap. Before we discuss the influ¬ 
ence of screening on the particle motion and radiation, 
let us consider phenomena in the gap. 

Aside from the initial moments when the particle is 
accelerated fast, it moves along equilibrium drift trajec¬ 
tory. In the equilibrium regime when the acceleration by 
electric field is balanced by energy losses of radiation, the 
position of the maximum of the energy distribution e* is 
given by (see Appendix [C]) : 

«.-(!)*(9i) 

where rg is the curvature of the magnetic field line, and e 
is the elementary charge. The equation is valid for a < 1. 
Note that in the derivation of Equation m, the energy 
losses were taken in the form given by Equation (l48l) . 
For comparison, if the cooling of electrons is dominated 
by synchrotron losses, we obtain the classical limit for 
maximum energy of radiation: 


y 4 a f ’ v ’ 

where a^ = e 2 /he « 1/137 is the fine-structure constant. 
The comparison of Equations m and m shows that 
the former depends on the strength and the curvature of 
the magnetic field, whereas in Equation (1921) the mag¬ 
netic field only enters through the ratio a cos 9 = 

In the case of a strong electric field (a = 1, 6 = 0°) 
and the dominance of energy losses by synchrotron radi¬ 
ation, the maximum energy of radiation is restricted by 
€ syn ~ 160 MeV. If the energy losses of electrons occur in 
the curvature or synchro-curvature regimes, the position 
of the maximum moves to higher energies. For example, 
for parameters used in our numerical calculations, the 
radiation maximum can be as high as 30 GeV (see Fig¬ 
ure [a. Thus, less intense energy losses of the curvature 
radiation lead to more energetic radiation. Remarkably, 
with an increase of curvature radius ro the radiation be¬ 
comes more energetic. 

For comparison of the drift approximation and the 
exact calculation, it is convenient, as above, to intro¬ 
duce the g-parameter. However, in cases when only the 
magnetic held presents, the g-parameter has been deter¬ 
mined by Equation (1561) as ratio of the curvature of the 
real trajectory and the curvature of the magnetic held 
lines. Now, when we compare the real trajectory and 
the electric drift trajectory, we determine g-parameter 
to be g = Re/Rc , where Re and R c are the curvature 
radii of the electric drift and the real trajectory, respec¬ 
tively. Then the characteristic frequencies of the radia¬ 
tion in both cases are related as e c = qeE ■ One can see 
from the right panel of Figure fl3l that in the gap (green 
lines) the Lorentz factors are close and the g-parameter 
is practically unity. The small difference in the Lorentz 
factors produces a small difference in the radiated spec¬ 
tra presented in the left panel of Figure [13] by red and 
green lines. This difference could be explained by the fact 
that the trajectories in the drift approximation and the 
exact calculation are not identical, as can be seen from 
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Figure 12. Projection on the XZ plane of t he trajectories of the 
electrons emitting the spectra in Figure 1131 The initial Lorentz 
factor is 70 = 100. The particles begin to be accelerated from the 
null surface at x = 0ARi c . The distances of axes in units of the 
light cylinder. The border of the gap is at x = 0.7 Ri c . The tra¬ 
jectories outside the gap are presented for the case of exponential 
screening with the width W = 0.1 R[ c . 

Figure El and the particle experiences slightly different 
magnetic and electric fields. 

When the particle escapes the outer gap it experiences 
the decrease of the field-aligned electric field. We ex¬ 
plore how this drop in strength influences the radiation 
of the particle depending on the intensity (gradient) of 
the drop. The relevant Lorentz factors and g-parameter 
are shown in the right panel of Figure [13] In the case of 
the slow decrease of the field-aligned electric field with 
W = 10 ~ 1 Ri c (pink lines), the g-parameter is unity and 
the Lorentz factors for approximate and exact calcula¬ 
tions follow each other, gradually decreasing. In this 
case the radiation spectra are close. For a more quickly 
decreasing electric field (W = 10 _2 I?z c , blue lines) the 
g-parameter oscillates around unity, and there is a slight 
divergence between the Lorentz factors calculated in dif¬ 
ferent approaches (especially close to the gap border), 
but the resulting spectra do not differ significantly. A 
different situation occurs when the field-aligned electric 
field drops very fast with screening width W = 10 -3 i?z c 
(brown lines). The g-parameter increases quickly on sev¬ 
eral orders of magnitude. At the same time, the Lorentz 
factor in the exact calculation drops, whereas in the drift 
approximation it continues to change gradually. This 
produces different spectra calculated in two approaches. 
The exact calculation gives more energetic radiated pho¬ 
tons compared to the drift approximation (solid lines in 
the left panel of Figure [0. Thus, in the case of fast 
screening, the drift approximation does not work. 

7. DISCUSSION 

In curved magnetic fields the particle trajectories ac¬ 
quire a complex structure when the curvature of the av¬ 
eraged drift trajectory is superimposed on the regular 
curvature of its gyration. The interplay of the two com¬ 
ponents produces a variable curvature of the particle tra¬ 
jectory. This results in the variability of radiation where 
the intensity and the spectrum are changed with the pe¬ 
riod of the particle gyration. The variability becomes 


pronounced during the transition from the synchrotron 
to the curvature radiation regime when the two com¬ 
ponents of the curvature of the trajectory become com¬ 
parable. During this transition the radiation cannot be 
described by synchrotron or curvature modes, but rather 
by their combination ■ synchro-curvature radiation. 

Because of the fast variability, only the averaged spec¬ 
trum of radiation can be observed. The spectrum aver¬ 
aged over the period of gyration is described by Equa¬ 
tion and depends on two parameters. At small or 
large values of the ?j parameter, the spectrum coincides 
with the standard form of the synchrotron spectrum. It 
is important to note that the radiation could not be re¬ 
duced to the formalism of the synchrotron radiation by 
introducing a single effective curvature radius. Indeed, 
the intensity of radiation averaged over the period of gy¬ 
ration is given by Equation (1501) . It can be represented 
in the standard form of Equation (1151) . if we introduce 
an effective curvature radius R e s = ro / \J 1 + r ] 2 . How¬ 
ever, if we want to have a standard exponential term e-x 
in an explicit form (as in the case of the synchrotron 
or curvature radiation) in the asymptotic presentation 
of the spectrum given by Equation (l64l) . the effective 
curvature radius should be defined differently, namely 
R e s = r 0 /{ 1 + 77 ). 

The conditions for the realization of one or another 
regime are determined by the relation between the drift 
velocity (3d and the velocity perpendicular to the drift 
trajectory /3±. The general conditions of the applicabil¬ 
ity of synchro-curvature radiation are also described in 
terms of Bn and /3± by Equation (l55l) . In the model case 
of a magnetic field of constant curvature and without 
energy losses, (3± is a free parameter. If (3± = 0, the par¬ 
ticle can move along a drift trajectory without gyration. 
This case is associated with the pure curvature radiation, 
which occurs solely due to curvature of magnetic field 
line. In the strong magnetic field, the drift trajectory is 
close to the field lines, so the substitution of the curva¬ 
ture of the field lines instead of the curvature of the drift 
trajectory is an acceptable approximation. In particular, 
in the radiation formulas we have neglected this differ¬ 
ence assuming that the condition given by Equation (0 
is fulfilled. 

The spectrum of radiation becomes quite sensitive to 
the pitch angles if (3± > /3d , with a shift of the maxi¬ 
mum to higher energies roughly as 1 + 77 . Therefore any 
angular distribution of particles wider than (3d produces 
a superposition spectrum with a high-energy cutoff that 
falls down slower than the standard exponential cutoff 
of the synchrotron or curvature radiation. Taking into 
account that in a strong magnetic field (3d is very small, 
even in the case of a very narrow angular distribution, 
the radiation spectrum could differ significantly from the 
unidirectional beam. This is demonstrated in Figures [2] 
and [3j where the radiation spectra is averaged over the 
pitch-angles of particles, assuming that Gaussian type 
and uniform angular distributions of particles are pre¬ 
sented. For the Gaussian type angular distribution, the 
spectral flux density can be presented in a simple analyt¬ 
ical form by Equation 0 , which provides good accuracy 
at x > 0.5. As seen from Figure [3l spectra becomes 
harder compared to the curvature radiation (£ = 0). It 
could also be noted that different angular distributions 
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Figure 13. Left panel: cumulative radiation spectra of electrons accelerated in the gap (red and green lines), and radiation spectra 
of electrons outside the gap (blue and purple lines) with screening width W = 10 ~ 1 Ri c (dotted lines), W = 10 ~ 2 Ri c (dashed lines), 
W = 10 SRic (solid lines) obtained in the exact numerical calculations (blue lines) and using the drift approximation (purple lines). Right 
panel: evolution of the Lorentz factor and (/-parameter with time. 


with the equal width give similar spectra. 

This interesting feature may be a key for the inter¬ 
pretation of the recent observations of pulsars by Fermi 
LAT that indicate that the energy spectra o f some pul- 
sars, in particular the Crab pulsar dBuehler et al.ll2012 ). 
in the cutoff region are significantly harder than e~ x , as 
predicted by the curvature radiation models. The energy 
spectrum of the puls ed emission of the Crab reported in 
(jBuehler et al.112012 ) can be readily described by Equa¬ 
tion (ITBll assuming £ = 1 and a quite reasonable value 
for the ratio 7 3 /ro — 10 13 cm” 1 . It should be noted that 
the form of the high-energy cutoff does not depend on en¬ 
ergy losses because it is produced by the most energetic 
particles at the first moments of radiation before they 
lose their energy. The angular distribution with enough 
broad angular width could appear when the particles ac¬ 
celerated along different field lines enter the region of the 
screened electric field. 

The results of the simplified analytical approach are 
confirmed by detailed numerical calculations. More im¬ 
portantly, the numerical calculations allow a rather gen¬ 
eral treatment of the problem extending to the case 
when the energy losses of electrons are taken into ac¬ 
count. We studied the radiation for the parameters typ¬ 
ical for the polar cap and outer gap models of the pul¬ 
sar magnetosphere following the evolution of the radia¬ 
tion of electrons with different initial pitch angles in a 
dipole magnetic field. As expected, the least energetic 
radiation is produced by electrons with an initial an¬ 
gle along the drift trajectory. The particle radiates in 
the synchro-curvature regime if the initial pitch angle is 
of the order of /3p, in particular, if the initial direction 
is along the magnetic field line. The spectra are sub¬ 
stantially more energetic when the initial pitch angle ex¬ 
ceeds (3d- For typical parameters of the polar cap model 
Pd = 2.2 • 10 ” 9 ( 7 / 10 8 ), so even a slight deflection from 
the field line results in the production of synchrotron ra¬ 
diation in the quantum regime. Moreover, this creates 
the interesting feature of a double peaked spectrum (see 
Figure 0) if the initial pitch angle is l/qo- 

In the context of the quantum synchrotron radiation, 


we should note the following circumstance. The pho¬ 
tons are produced in a quantum regime when the pitch 
angle of the charged particle is sufficiently large. The di¬ 
rection of the produced photons is basically along the 
direction of the charged particle. Thus, the photons 
are produced at large angles to the magnetic field. The 
fact that the photons have approximately the energy of 
the parent charged particles implies that the parameter 
X = -B 7 sin aj B cr > 1, thus the photon should produce 
an electron-positron pair and initiate an electromagnetic 
cascade. In this case, the cascade starts earlier than 
in the case of gamma-rays produced along the magnetic 
field lines. 

Due to the variation of the strength and the curva¬ 
ture of the dipole magnetic field, the change of radiation 
regimes is a quite complex process as can be seen from 
the variation of the g-parameter in Figures HE The 
transition to the curvature regime occurs with different 
rates for the polar cap and the outer gap environment. 
Whereas in the polar cap the particle makes a transi¬ 
tion to the curvature regime very fast (see Figure hod, 
the particle in the outer gap can radiate in the synchro¬ 
curvature regime most of the time while passing the gap 
(see the second panel of Figure 0. Although in the polar 
cap the synchrotron regime only lasts briefly, the particle 
loses most of its energy during this regime. Namely, if 
the initial pitch angle is a/70 (a > 1), the initial Lorentz 
factor on the drift trajectory is about 70/a. Since the 
typical Lorentz factor is 7 « 10 7 , the huge energy losses 
could be caused by a very small pitch angle. 

The treatment of pitch angles is not self-consistent 
without the consideration of particle acceleration. To 
study this question we examined the model of the accel¬ 
eration in the ele ctric a nd magnetic fields of the vacuum- 
retarded dipole (lDeutschlll955fl . as well as other simple 
models of the electric field. The calculations show that 
the particle basically moves along the drift trajectory 
jointly defined by the electric and magnetic fields. There¬ 
fore the appr oach of drift app roximation applied in the 
work of llvalapotharakos et all (| 2012 f) gives results that 
are similar to those obtained in the exact consideration. 
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The difference could arise when the particle escapes the 
acceleration gap experiencing the decrease of the field- 
aligned electric field. It has been shown that in the case 
of the fast screening of the electric field the particle can 
start to radiate more intensely. As a result, the parti¬ 


cle radiates more energetic photons than is predicted by 
the drift approximation. In general, the drift approxi¬ 
mation and the curvature radiation can only be applied 
to the equilibrium situations when the model parameters 
change slowly. 


APPENDIX 

MOTION OF A CHARGED PARTICLE IN ELECTRIC AND MAGNETIC FIELDS 

The motion of an ult rarelativistic cha rged particle in the electric (£) and magnetic (B) fields is described by the 
system of equations (lLandau fo Lifshit3119751 ): 


dr 

dt 

d(3 

dt 


= cP, 

e 


(Al) 


mcy 


(£-(/3£)/3 + /3xB), 


§ = A (W _ |fl 

dt me ime Ri 


where the curvature radius is expressed as 

1 

R r 


dp 


cdt 


77 TC 2 7 


^(£ + pxBY^(pey. 


(A2) 


The motion of an ultrarelativistic charged particle in the drift approximation is determined by the following system of 
equations: 


dr 

di = cf> - 

at me 6mc Hi 


(AS) 


where 


£*B B 

P B 2 


and / is defined from the condition that \P\ = 1 


f = \ 1- 


ExB 

B 2 


Then the curvature of the trajectory is 


Rc 


= \K\ = \{PV)P\. 


(A4) 


(A5) 


(A6) 


Note that the differential equations that describe the evolution of Lorentz factor have the same form in both approaches. 
The difference is in the determination of the curvature radius. 


ENERGY LOSSES IN THE QUANTUM REGIME 

In a strong magnetic field, the ultrarelativistic electrons can radiate in the quantum regime, provided that 

X = —qsina > 1, 

JDcr 


(Bl) 


where B cr = 2 ™ e fr « 2.94 • 10 13 G. In the presence of electric field, B sin a = |(/3 x B) \ should be substituted by 
yj{£ + P x B) 2 — ( P £) 2 . The energy loss rate can be written in the form (|Baver et alJll973t iBerestetskii et al.ll 19821) 

dE 


dt 


e 2 m 2 c 3 

\/ 3 irh 2 


H(X), 


where 


and 


H{x)= H(r,x)dT , 


H{t,x) = X [(1 - t)F(x) +xt 2 K 2 / 3 {x)\ , x = T _ , 


(B2) 

(B3) 

(B4) 
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where K 2 / 3 (x) is the modified Bessel function of the order 2/3, F(x ) is the emissivity function of the synchrotron 
radiation (see Equation [5U), r = e/E, where e is the energy of the radiated photon, E is the energy of the radiating 
particle. 

For calculations it is convenient to express H{\) in Equation (IB2I) in a simple approximate analytical form. Using 
asymptotics of this function 


H(X) 

H( X ) 


87T"\/3 2 
27 X ’ 
32ti 


X <1, 


243 


2 2/3 r 


(B5) 


X 


we suggest the following approximation 

8nV3 


H( X ) 


X 


27 


1 + 


3 (2x) 2/3 ] 

'(§)/ 


1 + 


2/3 , X»l- 

0.52^x(l + 3^-3.2x) 

l + 0.3yx + 17x+HX 2 


(B6) 


The first part of Equation (IB6I) (before the sign x) gives right asymptotics at y ^ 1 and x 1 and provide an 
accuracy better than 10% for other values of y, whereas the last term in the brackets improves the accuracy down to 
0.1% for any y. 

The spectrum of radiation in the quantum regime is described by the function that can be written as 


F q (x,r) = (1 - t)F(x) + t 2 xK 2 / 3 (x), x = 


E 


1 - r e c ’ 


(B7) 


where e c = 3eh ^ s ‘ n a y 2 is the characteristic energy of the emitted photon. To use this function in Equation (1571) . 
the variable e c sliould be replaced by frio*. An analytical approxi mation o f this function can be obtained using the 
approximation for emissivity function of the synchrotron radiation (lAharonian et ahll2010f) 


F(x) « 2.15ai 1/3 (l + 3.06a;) 1 / 6 


1 + 0.884a; 2 / 3 + 0.471a; 4 / 3 
1 + 1.64a; 2 / 3 + 0.974a; 4 / 3 ' 


and 


xK 2/3 (x) « 1.075a; 1 / 3 (l + 3.72a;) 1 / 6 


1 +1.58a; 2 / 3 + 3.97a; 4 / 3 


(B8) 


(B9) 


1 + 1.53a; 2 / 3 + 4.25a; 4 / 3 

Both approximations provide an accuracy better than 0.2% for any value of the argument x. 

MAXIMUM OF THE CURVATURE RADIATION SPECTRUM 

The maximum energy of the electron determined from the competition between the acceleration and the energy 
losses, is found from the condition t aC c = Uoss ■> where 


tloss 


E 


and 


tacc \dE/dt\ c 


\dE / dt\i oss 
E E 


3 mcR 2 
2e 2 7 3 ’ 


mc'y 


cecr cos 9B ea cos 6B 


(Cl) 


(C2) 


Note that the acceleration rate is defined as \dE/dt\ acc = ceS cos 0, where S is electric field, a = S/B is the ratio of 
the electric and magnetic fields, and 9 is the angle between them. From ti oss = t acc we obtain Lorentz factor of the 
electron: 

' 3cr cos 9BR ?' 


7eq — 


2e 


(C3) 


Substituting this equation into the expression for the characteristic energy of radiation given by Equation (1531) we find 

7 ' 3 

a cos 9B 


e* = 


Chy/K. 


(C4) 


Finally, taking into account that the curvature of the drift trajectory is R c = ro(l + (3%), where /3e = sin a is the drift 
velocity due to the electrical drift, and a is the pitch angle given by Equation (l89ll . we obtain 


e * = > 2 


3 \ 4 / a cos 9B\ 4 


chyjr 0 (l + ct 2 sin 2 9 ), 


(C5) 
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where rg is the curvature of the magnetic field line. 
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